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ABSTRACT 

The sharp magnetic discontinuities which naturally appear in solar magnetic 
flux tubes driven by turbulent photospheric motions are associated with intense 
currents. Parker (1983) proposed that these currents can become unstable to a 
variety of microscopic processes, with the net result of dramatically enhanced 
resistivity and heating (nanoflares) . The electric fields associated with such "hot 
spots" are also expected to enhance particle acceleration. We test this hypothesis 
by exact relativistic orbit simulations in strong random phase magnetohydrody- 
namic (MHD) turbulence which is forming localized super-Dreicer Ohm electric 
fields {Eu/Ed = 10^ ... 10^) occurring in 2. .15 % of the volume. It is found that 
these fields indeed yield a large amplification of acceleration of electrons and ions, 
and can effectively overcome the injection problem. We suggest in this article 
that nanoflare heating will be associated with sporadic particle acceleration. 

Subject headings: acceleration of particles — turbulence 

Understanding the mechanisms behind the dissipation of magnetic energy in the solar 
atmosphere is a key ingredient for the solution of several problems related to coronal heating, 
flares, and coronal mass ejections. Until recently, the study of magnetic energy dissipation 
seemed to follow two very distinctive paths: (1) phenomena related to the "Quiet Sun" 
and coronal heating have been interpreted as continuous wave dissipation (Hollweg 1984; 
Ulmschneider, Rosner & Priest 1991), (2) flares on the other hand were associated with 
"impulsive " dissipation. Most flow charts proposed for the second process start with the 
formation of current sheets in comparably simple topologies (e.g.. Priest & Forbes 2002), 
which reconnect, eject jets, and thereby drive turbulence. The turbulence, in turn, acts as a 
particle accelerator (Miller et al. 1997; Benz and Saint-Hilaire 2003), and finally dissipates 
into heat. 

Parker (1983) questioned the split of magnetic dissipation in waves (for the heating) 
and current sheet formation (for the flare). He proposed instead that random photospheric 
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footpoint motion forces the magnetic flux tubes to develop many tangential discontinu- 
ities throughout the corona, and pointed out that the associated currents, when exceeding 
a critical value, will drive local instabilities which rapidly release the magnetic energy in 
what he called nanoflares. Macroscopically, the instabilities manifest as localized anomalous 
resistivity. The work of Parker was followed by many articles analyzing how the photo- 
spheric motions couple to the corona (Heyvaerts and Priest 1984; Cargill 1993; Gudiksen 
and Nordlund 2002). On a separate development, numerical studies of decaying resistive 
MHD turbulence reveal the formation of intense localized current sheets (Matthaeus and 
Lamkin 1986; Biskamp and Miiller 2000), and simulations of nonhnear Alfven waves in a 
single magnetic loop (Moriyashu et al. 2004) show the sporadic occurrence of slow and fast 
shocks-mode shocks as dissipative discontinuities. Parkers ideas and the numerical studies 
seem to share one important aspect: the intermittency of localized currents inside the large 
scale structures. 

Most of the literature triggered by Parkers well-known conjecture has focused on the 
role of nanoflares in coronal heating. The electric flelds associated with the anomalous 
resistivity are, however, efficient particle accelerators as well. A simple yet realistic model 
for this kind of accelerator is homogeneous evolved MHD turbulence hosting intense localized 
current sheets. Models for particle acceleration using this scenario have been developed in 
past (Matthaeus and Lamkin 1986; Ambrosiano et al. 1988), and a recent article studies the 
non relativistic test particle motion in the electromagnetic environment of fully developed 
isotropic turbulence (Dimitruk et al. 2003) with uniform resistivity. 

Several observations seem to support the connection of heating with particle accelera- 
tion: the classical Lin et al. (1984) balloon observation of hard X Ray (HXR) microflares 
in active regions; ultraviolet (UV) subflares with HXR microflare counterparts (Porter et al. 
1995); tiny flares at centimeter wavelengths which are associated with Soft X Ray (SXR) 
transients (Gary et al. 1997); high-sensitivity observations of small decimetric reversed-type 
III bursts (Benz et al. 2001) suggesting downward electron beams and high-located accel- 
eration sites; acceleration without flares (Trottct 1994); SXR micro-events with associated 
gyrosynchrotron radiation from the quiet sun (Kruckcr et al. 1997); ubiquitous nano-events 
from the quiet sun observed in coronal extreme UV and radio radiation (Krucker and Benz 
1998); nonthermal tails in very small X-ray bursts (Krucker et al 2002; Benz and Grigis 
2002). In these observations, HXR and radio is believed to be a direct signature of nonther- 
mal electrons, while SXR and (extreme) UV are secondary effects after thermalization. 

In the present article we analyze the efficiency of particle acceleration in MHD turbulence 
with anomalous resistivity as a proxy for the solar corona. Unlike previous studies (Dimitruk 
et al. 2003; Gudiksen and Nordlund 2002) we assume that the resistivity is enhanced locally 
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at the places were the "hot spots" appear. 

Formation of hot spots. We consider coUisionless test particles in evolved homogeneous 
MHD turbulence with electromagnetic fields 

B = V X A (1) 
E = -d,A + ri{i)j, (2) 

where /^oj = V x B and r/Q) = 77o6'(|j| — jc) is an anomalous resistivity switched on above 
the critical current jc ~ enCg (Papadopoulos 1979). Here (n) the sound speed (number 
density) of the background plasma. The vector potential A(x, i) is modeled as a random 
field, subject to the MHD constraints 

E-B = if ri{j) = Q and E/Br^VA- (3) 

Equation (3) can be fulfilled in several ways. We use here a spectral representation in 
axial gauge, A(x, t) = ^^a(k)cos(k • x — uj{\<i)t — 0k) with a(k) • = and dispersion 
relation uj{\i) = k ■ v^, which is an exact solution of the induction equation with a constant 
velocity field v^. For simplicity, A(x, t) is taken as Gaussian with random phases 0k and 
(independent) Gaussian amplitudes a(k) with zero mean and variance 

(|a(k)|^)cx (l + k^Sk)-^ (4) 

A constant magnetic field Bq along can be included without violating eq. (3). The total 
MHD wave velocity is v\ = B'^{^qp)~^ with B"^ = Bl + (t| and (t| = | Xlk 1^ ^ ^(k)P the 
magnetic fluctuations. The matrix S = diag (/^, /^) in eq. (4) contains the outer turbulence 
scales, and the index v determines the regularity of the two-point function at short distance. 
The presented simulations have v — 1.5, = (0,0,va), and one turbulence scale is by 
an order of magnitude longer than the others, which describes migrating and reconnecting 
twisted fiux tubes (see Fig. 1 for an illustration). 

The vector potential contains some hundred wave vectors in the inertial shell min(/j"^)<|k|< 
10-2 -r^^ with the rms thermal ion Larmor radius. We focus on strong turbulence {gb/Bq 
> 1). The rms magnetic field B is a free parameter, which defines the scales of the particle 
orbits. The localized enhancement of the resistivity will (1) enhance the local heating inside 

the unstable current layer, Qj = rjjj'^, forming what we call here "hot spots". The fast 
heat transport away from the hot plasma will soon transform them to hot loops and (2) will 
dramatically enhance the particle (ion and electron) acceleration. The role of hot spots on 
coronal heating, their filling factor and their statistical characteristics will be analyzed in a 
separate publication. We focus in this article on the role of hot spots as particle accelerators. 
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The physical units used in this study are selected to represent the solar atmosphere. In 
SI units and for typical values B ~ 10~^ T, n ~ 10^^ m~^, T ~ 10^ K, the reference scales 
are as follows (electron values in brackets): time fi^^ ~ 10"^ s (6 ■ 10"^'^ s); length cQ~^ ~ 
300 m (0.17 m); thermal velocity ~ 1.2 ■ 10^ ms^^ (5 ■ 10^ ms^^); sound speed ~ 1 ■ 10^ 
ms~^; Alfven speed f a ~ 2 ■ 10^ ms~^; electron-ion collision time r ~ 0.003 s; Dreicer field 
Ed = ne^ In A/ {4:71 e^kTe) ~ 3-10~^ Vm~^. Time is measured in units offl~^=m/qB; velocity 
in units of the speed of light; distance in units of cQ~^. 

Particle dynamics. Particle momentum is measured in units of mc; vector potential in 
units of mc/q; magnetic field in units of B; electric current density in units of QB/{fioc)j so 
that the dimensionless threshold current is = {m/mp)csc/v\. The electric field is measured 
in units of cB, so that the dimensionless Dreicer field is (Vg/r')(me/m) with v'^ the electron 
thermal velocity and r' the electron-ion collision time. The dimensionless (primed) equations 
of motion are 

% = ^ (^) 

^ = v'xB'-f .VdJ'DJ' (6) 

with 7 the Lorcntz factor, B' = V' x A', and j' = V x B'. The dimensionless resistivity 
Tj' is characterized by the resulting Ohm field Eq^ = r;o|j| relative to the Dreicer field E^. 
Equations (5) and (6) are integrated numerically. 

When an initially maxwellian population is injected into the turbulent electromagnetic 
fields (eqns. 1-2), the particles can become stochastically accelerated. Figure 2a shows the 
energy evolution of protons with initial temperature T = 10^ K for the case rj' = 1.6 ■ 10~^, 
corresponding to Eq/Ed ~ 10^, which occurs in about 10% of the volume. The outer 
turbulence scales are = = 3 and l'^ — 60, and the rms magnetic field is B — O.OIT, 
with a background contribution Bq — O.OOIT along z. The density and temperature of the 
background plasma is lO^^m-^^ and T = lO^K, so that = 0.0004 and v'y^ = 0.007. For 
rj'—O and since a;<^Q, the particle motion is approximately adiabatic (injection problem). 
Finite resistivity 7]'>0 breaks adiabaticity, and energy can grow. The equation of motion 
(6) adds up - potentially - independent increments of kinetic momentum P' = 7v', so that 
this quantity may be expected to behave diffusive. In fact, after a short initial phase, (P'^) 
increases linearly with time (Fig. 2a, {P''^) ~ t). The momentum diffusion coefficient D 
(defined by (P'^) = A + Dt") increases, however, slower with rj' than Var{EQ) oc 1]'"^, and 
above En/ Ed - 5-10=^ subdiffusive behavior {a < 1) is observed (Fig. 2b, {P'^) ~ t°-^^). The 
average energy at fixed time increases then hke r)'^-^ (Fig. 3). Acceleration is dramatically 
increased for ions, which can reach GeV in less than 60ms. Standard diffusion (a— >1) is 
reached in the limit rj'^O. 
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Due to their large inertia, protons gain energy in relatively small portions. This is not 
so for electrons. The momentum evolution of coUisionless electrons of the high-energy tail of 
a maxwellian is shown in Figure 4. The plasma parameters are similar as in the proton case, 
but the maximum wave vector is somewhat smaller so that the volume fraction with | j | > jc 
is 0.07 only. Since electrons have much smaller Larmor radius, they follow the field lines 
perfectly adiabatically and gain energy only when dissipation regions are encountered (Fig. 
1). The orbits then exhibit large energy jumps (Fig. 4), so that a Fokker-Planck description 
is inappropriate. 

In order to gain physical insight into the interplay of electric acceleration and magnetic 
confinement we consider the space of the two approximate invariants of our model, energy 
£ and canonical momentum along the adiabatic direction, as illustrated in Fig. 5 (left). 
Here, y is the slowly varying direction, and the background magnetic field is represented 
by Ayo — Bqx. Color code represents a theoretical estimate on d{£')Q/dt — {rj' j' ■ v'), 
obtained by assuming constant particle density on surfaces of constant S and Py. The 
tremulous line in Fig. 5 (left) is a sample trajectory, whose energy evolution is shown in 
Fig. 5 (right). Red (blue) color indicates positive (negative) theoretical d{S')Q/dt. As can 
be seen, the theoretical (ensemble) estimate is statistically sharp enough to reflect in an 
individual trajectory. The underlying mechanism is purely geometrical: conservation of Py 
and £' restricts Ay to a band {Ay — PyY < — 1, and since Ay is positively correlated 
to Ey ~ — AAj^, the instantaneous value of d£li/dt' ~ E'yO'y can be guessed from Py and £' 
(Arzner et al. 2002). As a result, particles drift towards higher energy in the red domain of 
Fig. 5 (left) until they are scattered into the blue domain, where they sink to the boundary 
of the red domain. In spatially homogeneous turbulence this life cycle repeats indefinitely. 

We have investigated the effect of resisitive hot spots on coronal stochastic acceleration, 
with evolved MHD "turbulence" modeled by Gaussian fields. The hot spots form sporadically 
when the electric current exceeds a critical threshold. They do not only create localized 
Ohmic heating (as a possible coronal heating mechanism), but are also efficient particle 
accelerators. 

• For vanishing resistivity the ions are slowly accelerated (second order Fermi) and the 
electrons remain adiabatic. As resistivity increases at the hot spots, ions and electrons 
arc accelerated efficiently. 

• Acceleration is not sensitive to the type of low frequency MHD waves used as long 
as these are able to form tangential discontinuities and drive locally unstable currents 
which will enhance the resistivity. 

• The acceleration mechanism proposed overcomes the injection problem. 
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• Heating and acceleration may have a common origin. 

• Numerous well known observations can possibly be explained i.e. long lasting accel- 
eration, type III bursts before the flare or without flares, nonthermal X ray emission 
from microflares, nonthermal emission from the quiet sun. 

Several important aspects are still missing from our current analysis. The first is the 
non-Gaussian nature of real MHD turbulence, with formation of larger scale intermittent 
structures. Second, the local reconstruction of magnetic fields (Parker 1993) will drive non- 

lincarly stable discontinuities and create avalanches as it was proposed by Lu & Hamilton 
(1991) and Vlahos et al. (1995). Finally, the neglection of collisions of energetic particles 
breaks down on time scales of seconds in the solar corona. 

We suggest that the hot spots developed naturally inside a turbulent plasma can be any 
type of flare i.e. nanoflare, microflare or a regular flare, depending on the sharpness of the 
discontinuity and the value of the resistivity associated with the unstable current. Based on 
the above analysis we believe that reconncction is hosted by turbulence and not the opposite. 
Although we have envisaged parameters representing the solar corona, the mechanism may 
have applications to other astrophysical situations like turbulent jets, upstream of shocks, 
or turbulence in accretion disks. 

We thank A. Benz for helpful discussions. This work has received partial support from 
the European Research Training Network under contract No. NPRN-eT-2001-00310. 
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Fig. 1. — Localized dissipation regions |j| > jc (yellow), magnetic field lines (dotted), 
and electron sample trajectory (blue-red, encodes -Ekin) for = ly = 1 km, = 20 km, 
Brms=(2,2,10)-10-3T, Bq = 2-10-^T, En/E^~10^ and dissipative volume fraction 7%. The 
trajectory covers 2.6-10^ fit, and energetization occurs at dissipation regions. 
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Fig. 3. — Dependence of the average proton energy at fixed time on tfie magnitude of the 
Ohmic field. 
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Fig. 4. — Evolution of electron kinetic momentum for /j.=/j/=lkm and /2=10km. Initial 
velocities are from the tail v > 3 Vth of a maxwellian of 10^ K. Top: full population (blue), 
sample trajectory (red), and ensemble average (black). Bottom: Ohm field along the sample 
trajectory. 
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Fig. 5. — Left: acceleration and deceleration domains in {S', P^)-space. Color code represents 
a theoretical estimate of the ensemble-averaged energy drift d{£')n/dt. Right: simulated 
energy evolution. Red (blue) color indicates positive (negative) theoretical d{£')n/dt. 



